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Abstract 

, An improved mass-spring model for a dripping faucet is presented. The model is constructed based on 

\^ • the numerical results which we recently obtained from fluid dynamical calculations. Both the fluid dynamical 

I calculations and the present mass-spring model exhibit a variety of complex behavior including transition to 

y—i ■ chaos in good agreement with experiments. Further, the mass-spring model reveals fundamental dynamics 

^ ' inherent in the dripping faucet system. 
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O ■ §1 Introduction 
0^ ! 

' Motion of dripping water from a faucet is a well-known example of chaotic dynamical systems easily seen in 
^ I daily life. T^^^/^jl .V flYi ilVilVi'i n^ O'^r i'^ 'i'^i^l'^ h a.ve been performed to measure dripping time intervals between suc- 
cessive dropsMBQoBBLI'Qu'liSllil'OB A variety of dynamical behavior has been renoited, such as periodic 
, and nonperiodic attractors, period doubling, intcrmittcncy, boundary and interior crises,Lf Ej' Hopf bifurcation,B' 
O multiple stabilityB and hysteresis. LF Bifurcation is induced by a very small change of the flow rate, which is 
f-H a main control parameter of the system. Return map of the dripping time intervals shows one-dimensiijual 
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structure like tangled strings, in other words, map function is multi-valued for flow rates in some rangesoBE. 
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In addition, "unit structure" lifP which repeatedly ajipears as the flow rate is varied has been observed in a 
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bifurcation diagram in a wide range of the flow rate 

Although experimental results have been accumulated, no unified explanation exists for the complex phe- 
' nomena observed in the dripping faucet system. So far, most of theoretical analyses have been made based on 
simple dynamical models that are composed of a variable mass and a spring (mass-spring model) . The original 
mass-spring model which was proposed by Shaw0^ is described by 



dz\ , dz 



m- 



kz-^— + mg, (1.1) 



dt \ dt J dt 

dm „ 

— = Q = const., (1.2) 
where z is the position of the forming drop, m is its mass, and g,k and 7 are constant parameters. The mass 



increases linearly with time by eq. (1.2). When the position z reaches a critical point, a part Am of the total 
mass breaks away in proportion to the velocity at that moment. Although the model, being solved using an 
analog computer, simulated some experimental attractors, the parameter values (g, fc, 7), the flow rate Q and 
conditions at the breakup momertLoLa drop were not presented explicitly. 

Sanches-Oritiz and Salas-Britoll3'llZ|) numerically investigated the mass-spring model in some regions of the 
parameter space (Q, g, h), where h is another constant characterizing the falling mass as Am = hmz. 

Without solving eqs. (|]|) and (|r|) numerically, Austiru^ assumed an approximate equation of a damped- 
oscillation form: 

z(t) = At^/^ + B sin (iot^/^ + e'^\ (1.3) 
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He explained repeats of unit structure composed of period-one and period-two motion in a bifurcation diagram 
in terms of a feed back loop in the phase angle Qi. 

An approximate equation of another damped-oscillation form was used by d'Innocenzo and Renna.EJ-* They 
obtained bifurcation diagrams including the period-doubling cascade to chaos from both the approximate equa- 
tion and the numerical integration of eqs. (1.1) and (1.2). 

Although the mass-spring models exhibit complex behavior including a transition to chaos, several essential 
issues seem to be unclear: Are assumptions included in the models realistic? How to choose the parameter 
values? Can the models basically explain the dynamical features of real dripping faucet systems in a wide range 
of flow rates? 

Recently, Fuchikami et al. proposed an algorithm of fluid dynamical calculations for a dripping faucet. 
The algorithm is based on Lagrangian description instead of Eulerian one. A drop, which is assumed to be 
axisymmetric, is decomposed into many sliced disks and the dynamics of the fluid is described in terms of 
time evolution equations for disks under the influence of gravity, surface tension and viscosity. The numerical 
simulation succeeded in reproducing well not only the shape of drops, but also the complex dynamical behavior 
observed in experiments, including the unit structure of the bifurcation diagram. Moreover, the fluid dynamical 
calculations yield a lot of information on, for example, the spring constant, the breakup conditions etc., that is 
essential if one aims to construct a simplified model of mass-spring type. 

In this paper, we introduce another mass-spring model, an improved one reconstructed from the knowledge 
of the fluid dynamical calculations. Our mass-spring model can systematically explain various aspects of the 
complex dynamical behavior of the dripping faucet system. It will be clarified how the repeats of the unit 
structure appear in bifurcation diagrams. It will also be explained how a multi- valued map of the dripping time 
intervals is obtained. 

We first present results of our fluid dynamical calculations in §2, based on which the improved mass-spring 
model is constructed in §3. The last section will be devoted to a summary. 



§2 Results of the Fluid Dynamical Calculation 

We present here the results of our fluid dynamical calculations (FDC) of the dripping faucet system for a 
faucet of 5 mm in diameter. The flow rates are in the range 0.19 cm^/s< Q < 0.35 cm^/s, which correspond 
to drop rates in the range from 2.2 drops/s up to 6.0 drops/s. 
We have chosen the length, time, mass, and pressure units as 

lo = ^T/^, to={T/pgy/\ mo = pll, Po = . (2.1) 

For water at 20 °C, = 0.27 cm, to = 0.017 s, too = 0.020 g and Pq =270 dyn/cm^. The faucet radius was 
fixed at a = 0.916 (= 2.5 mm). The unit of viscosity was chosen as r]o = {pT^ / gY^^ . Then viscosity is ry = 0.002 
for water at 20 °C. 

A control parameter is the flow velocity kept constant at the end of the faucet (=at the top of the fluid), 
which is related to the flow rate Q as Q = t^o^vq. A breakup parameter cTcrit was introduced in the simulation. 
When the minimum cross-section of the necking region, CTmin = '"■'^min) reaches this value, we separated the drop 
from the remaining part and continued the simulation for the residue. (Another breakup parameter e deflned 
in the previous paper relates to CTcrit as e = CTcrit /'I'a^-) For saving tha computational time, a relatively large 
value was taken for the breakup parameter for long-term simulations namely, CTcrit was set to be 10~^ (i.e., 
eR:!4x 10-3). 

2. 1 Fluid motion of a dripping faucet 

To construct a simplified model, it is important to observe the fluid motion obtained from FDC carefully. 

Deformation of water drops is presented in Fig. |l| for the flow velocity vq = 0.07 and for the breakup parameter 
CTcrit = 10"^. Shape at a breakup moment is shown in Fig. |^(a). The breakup point of main drop appears at the 
lower root of a thin columnar bridge between the pendent mass (the conical part just below the faucet) and the 
nascent drop (the spherical part). Immediately after the separation of the main drop the thin columnar bridge 
recoils as the tip is formed like a round knob (Fig. |^(b)) and the secondary droplet, referred to as the satellite 
drop, is pinched off within a time At « 0.2 (Fig. [l](c)). The volume of a satellite drop is typically smaller than 
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Fig. 1. The fluid dynamical calculations (FDC). Time evolution of water drops for vq = 0.07. 



1 % of the largest drops. Just after the breakup moment, the fluid is stretched downward, so that the surface 
tension works as a restoring force. This causes oscillation of the fluid (Fig. [^(d)-(h)). When the increasing 
volume V of the fluid amounts to the static critical volume Vait (~ 4.7 for the faucet radius a — 0.916), thfi 
next process, "necking", begins, because the fluid with the volume larger than Vdit cannot be suspended.E^ 
In this process, a thin columnar bridge (between the pendant mass and the nascent drop) grows till the next 
breakup moment at which a round drop is separated from the remaining part. 

In summery, a breakup of a drop occurs through two main processes, the oscillating process and the necking 
one. Breakup of a small satellite drop might be neglected in a simplified model, because that occura in a 
relatively short time At w 0.2 after the breakup of a main drop. (As mentioned in the previous paperjij-* the 
size ratio of a satellite drop to the main drop becomes smaller for smaller faucet radii.) Dripping of a drop 
repeats through these two processes appearing alternatively. Figure |^ is a typical example for relatively slow 
flow velocities, i.e., for < vq < 0.16 when a = 0.916. In all numerical results appearing in the present paper, 
the above basic mechanism of breakup commonly works. 
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Fig. 2. Bifurcation diagram obtained from FDC. Plot of the dripping interval T„ vs. vqD 
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Fig. 3. FDC. (a) Plot of dripping time interval T„ vs. n for various values of vq. Value of vq is written in each 
frame, (b) Return map: plot of T„_|_i vs. T„. (c) Semi-log plot of power spectrum of (a) obtained from 2^ data 
points. 



2.2 Unit structure of bifurcation diagram 

Figure ^is a bifurcation diagram obtained from the simulation data of the dripping time interval T„, i.e., the 
time difference between the (n + l)-th and n-th breakup moments. We have ignored satellite drops that are 
smaller than 1 % of the largest drops. 

In Fig. ^, Period-one motion (PI), i.e., the motion in which the T„ value is unique, is seen repeatedly, for 
example, at vq = 0.0740 and vq — 0.0855. In a region between these two vq values, for example, 0.75 < vq < 
0.085, the Tn value is distributed over a finite range. This is seen in the form of a block, which repeats as vq 
increases. In the previous paper, we called this pattern unit structure. Looking at the bifurcation diagram with 
reference to T„, one sees each unit structure occ urs periodically, namely, at almost the same interval in r„. A 
similar periodicity is observed experimentally £3' 

Plot of Tn vs. n, its return map (i.e., plot of Tn+i vs. T„) and power spectrum are given in Fig. ^for several 
values of vq belonging to the same block, i.e., the same unit, whose parameter region is denoted as R5. 

Period-one motion at vq = 0.0997 period-doubles backward to P2 motion at vq = 0.0957. The latter further 
period-doubles backward at vq — 0.0951, which can be seen as a peak at the frequency / = 1/4 in the power 
spectrum of T!„. These results, together with the results for vq = 0.0948 in Fig. |[ suggest a backward period- 
doubling route to chaos. In each unit in Fig. ^, backward period-doubling from PI to P2 (as vq decreases) 
appears. In other words, the backward bifurcation from PI to P2 occurs repeatedly with very small intervals 
of Wo, and in each interval a series of backward period-doubling cascade to chaos may occur. This suggests that 
experiments of measuring the dripping time intervals can be greatly influenced by a very small uncontrolled 
change of the flow rate. Also it might be possible that some experimental data of PI, P2 and P4 motion which 
arranged as if they were a part of a series of period-doubling cascade may belong to different series of cascade. 
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2.3 Relation between the unit structure and oscillation of fluid 



We will show that the repeating unit structure in the bifurcation diagram is closely related to the oscillation 
of fluid. After a breakup, the position zq, the center of gravity of fluid, oscillates like a damped oscillator. 
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Fig. 4. FDC. Oscillation of the center of gravity of the fluid. Period-1 motion. Arrows indicate the moments 
at which V = Krit = 4.7. 

Let us observe first PI motion in detail. Figure ^ shows the oscillations of zq for three kinds of PI motion for 
vo — 0.0740, 0.0855 and 0.0977. Corresponding to these vo values, the periods (= the dripping time intervals) 
are T — 21.3, 18.7 and 16.6, which belong to three different units, whose regions in vq direction are i?7 , _/i!q 
and i?5, respectively. As mentioned in subsection p.l| , a breakup of a drop occurs through two main processes, 
oscillation and necking, which can well be observed in Fig. ^. Discontinuous changes in zq seen in Fig. || occur 
at breakup moments, at which a part of the fluid is lost. Arrows in Fig. ^ indicate the points at which the total 
volume V reaches 4.7 — Vaix. that is the static critical value. Let N be the oscillation frequency of zq the 
number of peaks in the plot of zq vs. i) during a dripping time interval. In Fig. 0, one sees N is 7, 6, and 5, 
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corresponding to which the regions of units are R-j, Rq and i?5, respectively. In this way, it was found that each 
unit is characterized by an integer N such that the unit includes PI motion whose oscillation frequency during 
a dripping time interval is N. 

After the volume exceeds the critical value, the fluid cannot be suspended and zq drops rapidly because of 
necking. Time of the necking process, Ainock, is not negligible in the sense that the total volume still increases 
for this time even after it reaches VcHt- In Fig. ^, A^neck, is almost equal to or more than 20 % of the total 
dripping time interval. For vq — 0.0977, the total volume at the breakup moment amounts to Vtotai ~ 5.8, about 
27 % of which is left after breakup. The necking time Ai„cck can change depending on details of the oscillation 
of the fluid, which causes variation of dripping time intervals. As an example, let us discuss the following P2 
motion. 
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Fig. 5. FDC. Same as in Fig. |[ Period-2 motion. 

In Fig. H, P2 motion is plotted for vq — 0.0705 and vq — 0.0798, which belong to two different units in Fig. ^ 
specified by the intervals Rj and R^, respectively. For vq = 0.0705, oscillations with N — 7 (the dripping time 
interval being Ti = 21.9) and N = 8 (the dripping interval being T2 = 24.0) occur alternatively. The last (= the 
eighth) peak of zq appears in every second time interval and is not so sharp as other peaks. Figure |^ shows that 
if the velocity zq at the moment of ^ = 4.7 = Vent (indicated by arrows) is upward, an extra weak oscillation 
occurs before breakup. Observing the shapes of drops at various stages, we found that the waist of the fluid 
(at which r = r^nin) is more slender when zq at the moment of V = Vent is downward than when it is upward, 
which causes a longer time interval T2 than Ti. 
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In summary, each unit in the bifurcation diagram is specified by the symbol Ri indicating a region of vq. The 
sufhx i equals the minimum value of N, oscillation frequencies of the fluid during the time interval between two 
successive breakup moments. 




Fig. 6. FDC. Same as in Fig. |^ but for different values of vq. 

2.4 Multi-valued return map 

Intermittent P3 motion is observed at vq = 0.0890, the spectrum of r„ for this motion showing a peak at 
/ = 1/3 (Fig. ^). A small change of the vq value makes the spectrum broad, which indicates chaotic motion, 
for example, at vq = 0.0888 and vq = 0.0892. In Fig. ^(b) we see that return maps for these vq values apijear 
to be multi- valued, at least two- valued. This kind of multi- valued maps are often observed in experimental!' ES' 
as well as in simulations.tHllP 

Figures |^(a) and 0(e) show other multi-valued return maps for vq = 0.0886 and vq = 0.0697. So far, we have 
plotted return map for dripping time intervals r„, i.e., {Tn^i vs. T„). Note that experimental maps are usually 
plotted in the same way because measuring other quantities like mass is not easy in experiment. To investigate 
the mechanism of the multi-valued map, we plotted return maps using data of mass, where the mass is equal 
to the volume because p = 1- Figures 0(b) and 0(f) are plot for the dripping mass, i.e., {md,n+i vs. m^.n)- 
They are again multi-valued maps. In contrast, the map function for the remnant mass, i.e., (m^^n+i vs. mr,n), 
is roughly single valued, although that is double humped (Figs. 0(c) and 0(g)). Furthermore, if we plot the 
dripping time interval vs. the previous remnant mass, i.e., {Tn^i vs. rrir^n), the return map becomes a simple 
unimodal function as seen in Figs. 0(d) and 0(h). This implies that the time interval Tn+i is directly determined 
only by the previous remnant mass mr.„. 

Relation between the single-valued map for mp,™ and the two-valued map for r„ is understood as follows. Let 
us consider two variables x and ?/, and assume logistic map for x: 

Xn+l ^ F{Xn) ^ AXn{l ~ Xn). (2.2) 
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Fig. 7. Return maps obtained from FDC. The upper row for vq — 0.0886. The lower row for vq = 0.0697. 
(a),(e) Plot of r„+i vs. T„. (b),(f) Plot of md,n+i vs. md,n where 77id,n represents the mass of the n-th falling 
drop. (c),(g) Plot of nir.n+i vs. mr,n where TOr,n is the n-th remnant mass. (d),(h) Plot of T„+i vs. mr^„. 




Fig. 8. A two- valued map (a) is obtained from a single- valued map (b) via a nonlinear asymmetric transforma- 
tion (c) from Xn to yn+i- 

Further, assume that the next value of y is uniquely determined by the present value of x as 

3\/3 

yn+l = G{Xn) = —^(1 - Xn)Xn{l + X^) (2.3) 

for < a;„ < 1 and < 2/n < 1- Then the map from ?/„ to j/n+i is given by 

yn+i = G{xn) = G{F{xn-i)) = G{F {G-\yn))) , (2.4) 

which is two- valued as illustrated in Fig. ||(a). Note that the transformation G{x) is asymmetric with respect 
to a; = 1/2. This asymmetry causes the multi- valued map for ?/„. 

The situation is roughly the same for the dripping faucet system. That is, x„ and ?/„ correspond to mr.„ and 
r„, respectively, which leads to a multi- valued return map for T„. 

Some attractors of a dripping faucet look similar to Henon map as shown by Shaw.ld^ However, the above 
explanation suggests a difference in the fractal structure of attractor between the present system and Henon 
map. The difference arises because the present system is basically described in terms of a one-dimensional map, 
whereas Henon map is two-dimensional. We will discuss this point more detailedly in the next section. 
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Fig. 9. Hysteresis observed in FDC. Plot of the dripping time interval r„ for the n-th drop, (a) Increasing flow 
velocity, (b) Decreasing flow velocity. 

2. 5 Hysteresis 

In Fig. 1^ is presented hysteresis: Different behavior is observed for some vq values depending on whether the 
vo value has been increasing (Fig. ^(a)) or decreasing (Fig. ||(b)). Hysteresis is observed also in an experiment 
by Sartorelli et a/.U In the next section, we will analyze our improved mass-spring model and show that the 
range of vq on which the hysteresis occurs is nearby the PI motion belonging to the border of units. 



2. 6 Mass dependence of the spring constant 

We have already seen that the fluid oscillates until the volume reaches the critical value T4rit- This justifies 
the mass-spring model. However, for the model to be more realistic, one needs to know its spring constant. 

To estimate the spring constant k, we did a test simulation of FDC as follows: We started the simulation with 
some small initial volume (Vq — 1.4) and some slow flow velocity [vq — 0.07). When the mass {— the volume) 
reaches a certain value m for which we want to know the spring constant, we stopped the influx of water (i.e., 
was reset at zero). Then the fluid oscillates with very weak damping. Figure |l^(a) shows the oscillation of the 
center of gravity of the fluid for several values of m. Measuring the period T(m) of oscillation after the influx 
of water is stopped, the spring constant was obtained as fc(m) — m{2Tr/T{m))'^. The mass dependent spring 
constant thus obtained is plotted in Fig. 10(b), which becomes zero near a critical value m (« Vcrit = 4.7). 
Based on these data, we will introduce a linear function k{m) in our mass-spring model (eq. ( |3.4| ) in |§3| ). 

2. 7 Summary of the fluid dynamical calculations 

Several aspects observed in the dripping faucet system are summarized as follows: 

1. Breakup of a drop occurs through two processes, oscillation and necking. 

2. Although the return map Tn+i vs. r„ is a multi-valued function, the return maps mr,n+i vs. mr_„ and 
T„+i vs. rrij-^n are single-valued. 
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3. Repeating unit structure is observed in the bifurcation diagram over a wide range of the flow velocity. 
Each unit includes period-doubling and tangent bifurcation in a narrow range. 

4. Each unit is specified by the oscillation frequency of position zq during the time interval Tn- 

5. hysteresis. 



§3 Improved Mass-Spring Model 

We propose an improved mass-spring model based on the results of FDC for relatively low flow velocities 
vo < 0.16 (i.e., Q < 0.35 cm'^/s or equivalently, 6.0 drops/s). 

We have observed zq, i.e., the movement of the center of gravity of the fluid in detail. Our mass-spring model 
describes only this degree of freedom explicitly like other mass-spring models. It was found from the FDC that 
for low flow velocities breakup generally occurs through the process of "regular" necking. Here, "regular" means 
that the necking begins after but not before the total volume has reached the static critical value V^rit • We take 
into account this issue in our mass-spring model. 

After breakup, the fluid starts to oscillate because of the surface tension which acts as a restoring force. The 
position Zq oscillates like a damped oscillator. Its mass-dependent spring constant has been obtained from the 
FDC. The dripping time interval strongly depends on the direction of the velocity zq at the moment at which 
V = Vcrit- If Zq is upward, an extra oscillation of zq occurs before the next breakup, which makes the dripping 
time interval longer. 



3. 1 Equation of motion 
Our mass-spring model obeys the following equations of motion; 



dp dz 

- = -kz-^-+mg^ (3.1) 
dp d?z f dz \ dm 



= m 



dt dt^ \dt V dt 



■li-^o (3.2) 



where 7 and g are constant parameters. The units in eq. ( 2.1 



have been used, so that g = 1. The damping 



parameter 7 was chosen as 7 = 0.05. The second term in cq. (^^) represents adhering of mass with a relative 
velocity (i — vq). The mass increases with a constant flow rate Q: 

dfn 

= Q = const. (3-3) 

The term vom (= vqQ) in eq. ( |3.2| ) was neglected because this is proportional to a small quantity Vq. (Note 
that Q = ira'^vo.) 

Referring to the results of Fig. |l^(b) , we assumed the mass-dependent spring constant k takes the following 
linear form: 

, r -11.4m + 52.5 (77K4.61) 
^("^) = | (m>4.61) ■ (3-^) 

When the mass m amounts to 4.61, the spring constant becomes zero, then the mass undergoes a free fall, 
which corresponds to the necking process of the fluid. An important difference of the present model from other 
mass-spring models is that the necking process is included explicitly. 

We have made several assumptions for reset conditions at breakup moment. It is assumed that breakup 
occurs when the position z reaches a critical point Zdit- A part of the mass is then lost and the mass is left. 
The remnant mass was renewed as 

nil- — 0.2m + 0.3, when z = Zcritj (3-5) 



where m is the total mass at the breakup moment. The linear relation in eq. (3.5) is oversimplified because 
this yields a single- valued map of md,n, which contradicts the FDC data (Figs, ^h) and 0(f)). If eq. ( |3.5| ) is 
replaced by a nonlinear function of m, then the md.n map also becomes multi-valued as well as T„ map. We 
shall revisit this issue later on. 
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The position and velocity are renewed as 



z = zo = const, 
i = 0, 



when z = ZcHt- i'^-^) 



The above reset conditions for z and z are also simplified approximations. However, we found from the data of 
the FDC (see Figs. H a nd that these approximations are, at least qualitatively, acceptable. 

Equations ( |3.5| ) and ( |3.6| ) imply that a set of the points on which the state is reset is located on a straight 
line in the phase space (z, z, m), which will be shown in Fig. |l^ later. We chose Zcrit = 5.5 and zq — 2.0. Note 
that the reset state depends only on the total mass just before the breakup. This implies the dripping time 
interval Tn+i is decided uniquely, (by solving eqs. (|3.l|) and (|3.2|)) from the previous value mr,™ of the remnant 



mass. Thus TOi-.k+i is uniquely decided only by mi..„ from eq. (3.5) because the total mass m at the breakup 
moment is also decided uniquely by mr.n via T„+i. In this way, the present model is described in terms of a 
one-dimensional map from mi..„ to rrir^n+i, which is consistent with data obtained from the FDC (Figs. |^(c) 
and @(g)). Return maps of the present model will be presented later (Fig. 

3.2 Bifurcation diagram 



We numerically solved eqs. ( |3.1| ) and ( |3.2| ) together with eqs. ( p.3| ) ~ ( |3.6| ) and obtained a bifurcation diagram 
over a wide range of flow velocities vq (Fig. |ll|). We chose a = 0.916. The diagram was made by raising vq. 
The final state at the present vq value is chosen as the initial state ^LJhfi next vq value. The repeating unit 
structure in bifurcation diagrams which was observed in experimentsoE^'liJ) and our FDC is reproduced. 

In Fig. |l^, bifurcation diagrams are shown for three different ranges of vq included in Fig. [l^. In a range of 
very small vq (Fig. ^(a)l, p^rh unit is composed of only PI and P2 motion. A similar type of unit structure is 
observed experimentally Bt^llJ We have not obtained this type of unit structure from FDC because long-term 
simulations for such small flow velocities take too long computational time to obtain enough data at many vq 
values. 

For higher flow velocities (Fig. |l^(b)), each unit includes forward and backward period-doubling cascades 
that start from P2 motion, in addition to the normal (i.£._ starting from PI motion) period-doubling cascade 
to chaos. Similar data are also obtained in experiments.HQ'lid^ The numbers in Fig. |l^ denote the oscillation 
frequency N during the dripping time interval r„ for the PI motion at the points indicated with allows. The 
frequency N increases by one by moving to the neighboring unit on the left, i.e., the side of lower flow velocity. 

The oscillation of the position z is represented in Fig. [ij. Figures |T3|(a) and |l^(c) show PI motion at two 
Vq values, Vq = 0.092 and 0.110, having N = 5 and 6 respectively. They belong to two neighboring units in 
Fig. ^2|(c). In P2 motion at vq = 0.100 (see Figs, ^(c) and p^(b)), two different frequency N = 5 and 6 appear 
alternatively. These features have been observed in the FDC (Figs. || and 0). 

Further, the present model exhibits hysteresis as seen in Figs. |l^(c) and p2|(d). Figures [T^(c) and |3|(d) show 
different types of PI motion observed at the same vq value: vq = 0.110. They belong to the lower and upper 
branch of hysteresis curve and obtained by raising and reducing vq, respectively. In Fig. |l^, circles indicate 
the moment at which the mass reaches the critical value m — 4.61 and the spring constant k becomes zero. 
Comparing Figs. |l^(c) and |3|(d), we see that the dripping time interval r„ is longer in (d) mainly because the 
velocity z is upward at the moment of fc = 0. The hysteresis will be discussed in more detail later on (in the 
next subsection) together with map functions. 

The repeating unit structure of the bifurcation diagram shows a scaling law, which has been found experimen- 



tally .HB'EJ^ As was mentioned in subsection |2.2| , each unit occurs periodically, i.e., at almost the same interval 
in T„. Let us specify bifurcation points of forward period-doubling from PI to P2 as (wo.jv, Tm^at), where N is 
the oscillation frequency. Then these points depend on N as 

r„,jv ~ 2.2N + 3.9, (3.7) 
vo^N ~ 0.58N-^°. (3.8) 



Equation (3.7) corresponds to the periodicity of the repeating unit structure in T„ direction. The scaling: law, 
eq. ( |3.8[ ), agrees well with the experimental scaling law Q (x N^^-^ obtained by Katsuyama and Nagata.Ej^ An 
empirical scaling law which agrees with the above for sufficiently large N was also obtained by Shoji.O 
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3.3 Return map 

The present mass spring model is basically described by a one-dimensional map. Once mi-,„ (the n-th remnant 



mass just after breakup) is given, T„+i (the next dripping time interval) is decided uniquely by eqs. (3.1) and 
( |3.2| ). The map from m^^n to T„_|_i, T„_|_i — G(mr,„), is shown in Fig. p^(c). Then the next value of the remnant 
mass, rrii-^n+i, is uniquely determined as 

mr,„+i = 0.2 (m,,„ + QG{m,,n)) + 0.3, (3.9) 



because of eqs. (3.3) and ( |3.5| ). (Fi g. |l4| (b).) Therefore, the map from r„ to Tn+i is also one-dimensional even 
if it can be multi- valued. (See Fig. |14|(a).) We have illustrated how the map can be multi- valued in Fig. ^. A 
similar return map of Tn+i vs. r„ which looks like tangled strings was also observed in the FDC (Figs, ^a) and 
(e)). This suggests the possibility that experimentally observed return maps of r„ can also be one-dimensional, 
even when the maps look like complicatedly tangled strings. 

The present mass-spring model is thus basically different with Henon map in the structure of chaotic attractor. 
This is because Henon map is two-dimensional, i.e., the next value Xn+i is decided not only by the present value 
Xn but also by the previous value Xn-i- 



-1 



= l-a< + fox„_i. (3.10) 



In Fig. |T5| some attractors are represented for the MSM and Henon map. At first sight, they look similar. 
However, Henon attractor includes selfsimilar structure on a two-dimensional plane, whereas the attractor of 
MSM is restricted on a one-dimensional manifold. Note that the situation does not change even if the linear 
reset condition (eq. ( |3.5[ )) is replaced by a nonlinear relation. The m&-p function Tn^j ^j-i-i vs. vricj then becomes 
multi-valued as seen in Figs. 0(b) and ^(f) of FDC but the dynamics is still described by a one-dimensional 
map, namely the next value Xn+i is determined only by the present value Xn- 

We can also improve the reset condition z — such that the renewed i value depends on the z value just before 
the breakup. Then the dimensionality of the system increases and the chaotic attractor may exhibit selfsimilar 
structure as Henon map. Such a refinement of the model is, however, not so important in the sense that the 
essential aspect, i.e., the return map which looks like tangled strings can be obtained by a one-dimensional map. 

Figure ^ shows map function TOr,Ti+i = -M(™r,ra,i^o) for several values of the flow velocity vq. The map 
M is an oscillating function of mr.n- The oscillation causes the repeats of unit structure. The amplitude and 
the maximum slope of the map function becomes large with vq. This is the reason why the unit structure is 
simple for small vq values and becomes more complicated as vq increases. In a range of small enough vq values, 
each unit includes PI and P2 motion only. In a range of larger vq values, period-doubling bifurcation can 
occur more than once in the same unit. For large enough vq values, each unit includes forward and backward 
period-doubling cascade to chaos and tangent bifurcation. 

Hysteresis can be understood from Fig. |l^(b). There are two stable fixed points at, for instance, va = 0.111. 
If this Vq value has been reached by an increasing process of vq, then the fixed point on the left (with a smaller 
TOr,n value) is realized. As vq increases further, the fixed point on the left shifts to the right (see for vq — 0.112 
in Fig. |l^(b)) until this fixed point disappears by forward tangent bifurcation. At this moment, the other fixed 
point (with a larger mr^„ value) is still stable (at a vq value slightly larger than 0.112) and this fixed point is 
realized, which causes discontinuous increase in the mr_„ value. In an decreasing process of vq, the fixed point on 
the right (with a larger mr.„ value) is realized first, and when this fixed point disappears by backward tangent 
bifurcation, the other stable fixed point on the left is realized. The hysteresis occurs in this way. 

Some experimental data have been explained in terms of boundary crisis.™ However, as suggested from the 
bifurcation diagram, our mass-spring model yields no boundary crisis. Boundary crisis occurs for example in 
a quadratic map when an unstable fixed point that is a basin boundary of a chaotic attractor collides with an 
edge of this chaotic attractor. In contrast, the unstable fixed point expressed by the dotted line in the inset of 
Fig. |l2| (d) is a basin boundary of two fixed point attractors and it disappears together with one of these fixed 
point attractors by tangent bifurcation. 

Finally we present chaotic trajectories in the phase space (z, i, m) in Fig. |l^. The state point is on the plane 
z = z^ciit = 5.5 at each breakup moment. A set of state points just after the breakup moments is on a line 
segment satisfying z = 2 and i = (reset conditions). This line segment is stretched and folded on the plane 
of z = Zcrit- Then the image again mapped onto the initial line segment by the reset conditions expressed as 



eqs. ( |3.5| ) and ( |3.6D 
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§4 Summary 



We introduced an improved mass-spring model based on the results of our fluid dynamical calculations. The 

FDC data reproduced the repeating unit structiu'e in the bifurcation diagram, stable P2 motion included in the 
middle of each unit, forward and backward period-doubling starting from this P2 motion, hysteresis etc. which 
were observed experimentally. The FDC showed that the repeating unit structure is caused by the oscillation 
of fluid between two successive breakup moments. The FDC also showed that the necking process plays an 
important role in changing the dripping time interval. 

Our mass-spring model corresponds to a one-dimensional map system, i.e., the simplest system that can yield 
chaos. In spite of that, the model succeeded in reproducing systematically a variety of dynamical behavior 
observed in many experiments. Especially, the model supplemented a part of the bifurcation diagram in a 
range of very small flow velocities vq that was not obtained by FDC because of long computational time. For 
these small vq values, the model showed that each unit structure in the bifurcation diagram is composed only 
of PI and P2 motion. For relatively large vo values, each unit includes period-doubling cascade to chaos. The 
period-doubling cascade starts not only from PI motion but also from P2 in the middle of the imit. For the 
latter, the cascade occurs both forward and backward. These kinds of structure agree well with experimental 
observations. The two- valued map for the dripping time interval was explained from a single- valued map of the 
remnant mass. Mechanism of the hysteresis was understood based on the return map. 

The simple dynamics of the present mass-spring model clarifled fundamental mechanisms inherent in the 
complex behavior of the dripping faucet system. 
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Fig. 10. FDC. (a) Oscillation of the center of gravity for several values of mass smaller than the critical value, 
(b) Mass dependence of spring constant. 




Fig. 11. Bifurcation diagram obtained from our mass-spring model (MSM) for a wide range of vq. The diagram 
was made by raising vq- The final state for the present vq value is chosen as the initial state for the next vq 
value. 
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Fig. 12. MSM. Bifurcation diagrams. Enlargement of several parts of Fig. ^ The number in each frame 
represents the frequency A'' of oscillations between two successive breakup moments. In (a), (b) and (c), the 
control parameter vo was raised, (d) is the same as (c) but vo was decreased. For (c) and (d), the bifurcation 
diagrams are different in some ranges of vq (indicated with <->), which shows hysteresis. The inset shows a 
hysteresis curve, in which a dotted line represents an unstable orbit. 
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Fig. 13. Oscillation of the position z in MSM. Circles indicate the moments at which the spring constant 
becomes zero, (a) PI motion, (b) P2 motion. (c),(d) Different types of PI motion at the same vo value: 
vo = 0.110. 
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Fig. 14. MSM. uo = 0.115. (a) Map for the dripping time interval T„, (b) for the remnant mass mr,^. (c) Map 
from mr,„ to T„+i. 
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Fig. 15. (a) MSM. Return map of r„ for vq = 0.094. (b) Henon map: Xn+i = -0.33yn + 1 - 2.07a;^, j/„+i = Xr, 



1.5 



Vo = 0.110 



1.2 




1.2 



1.5 



v„ = 0.111^ 
1)0 = 0.112- 
Do = 0.113 — 



1.2 

1.5 1.2 




• Stable fixed point 
o Unstable fixed point 



1.5 



Fig. 16. MSM. Map function mr,n+i = M{mr^n,vo) for various values of vq- 
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Fig. 17. Phase space portrait of MSM. vq = 0.115. The inset is a Poicare surface of section at the critical 
position z = Zcrit = 5.5. 
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